Targeted amplicon sequencing of wastewater samples for detecting SARS-CoV-2 variants with high sensitivity and resolution.

Wastewater-based epidemiology (WBE) is a promising approach for monitoring the spread of SARS-CoV-2 within communities. Although qPCR-based WBE is powerful in that it allows quick and highly sensitive detection of this virus, it can provide limited information about which variants are responsible for the overall increase or decrease of this virus in sewage, and this hinders accurate risk assessments. To resolve this problem, we developed a next generation sequencing (NGS)-based method to determine the identity and composition of individual SARS-CoV-2 variants in wastewater samples. Combination and optimization of targeted amplicon-sequencing and nested PCR allowed detection of each variant with sensitivity comparable to that of qPCR. In addition, by targeting the receptor binding domain (RBD) of the S protein, which has mutations informative for variant classification, we could discriminate most variants of concern (VOC) and even sublineages of Omicron (BA.1, BA.2, BA.4/5, BA.2.75, BQ.1.1 and XBB.1). Focusing on a limited domain has a benefit of decreasing the sequencing reads. We applied this method to wastewater samples collected from a wastewater treatment plant in Kyoto city throughout 13 months (from January 2021 to February 2022) and successfully identified lineages of wild-type, alpha, delta, omicron BA.1 and BA.2 as well as their compositions in the samples. The transition of these variants was in good agreement with the epidemic situation reported in Kyoto city during that period based on clinical testing. These data indicate that our NGS-based method is useful for detecting and tracking emerging variants of SARS-CoV-2 in sewage samples. Coupled with the advantages of WBE, this method has the potential to serve as an efficient and low cost means for the community risk assessment of SARS-CoV-2 infection.

Ensitrelvir is effective against SARS-CoV-2 3CL protease mutants circulating globally.

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a public health concern worldwide. Ensitrelvir (S-217622) has been evaluated as an antiviral treatment for COVID-19, targeting SARS-CoV-2 3C-like protease (3CLpro). Ensitrelvir has been reported to have comparable antiviral activity against some of the SARS-CoV-2 variants: alpha, beta, gamma, delta, and omicron (BA.1.18). In this paper, we describe that ensitrelvir is effective against newly emerging SARS-CoV-2 variants and globally prevalent 3CLpro mutations. Ensitrelvir exhibited comparable antiviral activity against SARS-CoV-2 variants, including recently emerging ones: omicron (BA1.1, BA.2, BA.2.75, BA.4, BA.5, BQ.1.1, XBB.1, and XE), mu, lambda, and theta. Genetic surveillance of SARS-CoV-2 3CLpro, the target of ensitrelvir, was conducted using a public database and identified 11 major 3CLpro mutations circulating globally (G15S, T21I, T24I, K88R, L89F, K90R, P108S, P132H, A193V, H246Y, and A255V). The 3CLpro mutation from proline to histidine at amino acid position 132 was especially identified in the omicron variant, with prevalence of 99.69%. Enzyme kinetic assay revealed that these 3CLpro mutants have enzymatic activity comparable to that of the wild type (WT). Next, we assessed the inhibitory effect of ensitrelvir against mutated 3CLpro, with it showing inhibitory effects similar to that against the WT. These in vitro data suggest that ensitrelvir will be effective against currently circulating SARS-CoV-2 variants, including omicron variants and those carrying 3CLpro mutations, which emerging novel SARS-CoV-2 variants could carry.

Mechanistic roles of microRNAs in hepatocarcinogenesis: A study of thioacetamide with multiple doses and time-points of rats.

Environmental chemicals exposure is one of the primary factors for liver toxicity and hepatocarcinoma. Thioacetamide (TAA) is a well-known hepatotoxicant and could be a liver carcinogen in humans. The discovery of early and sensitive microRNA (miRNA) biomarkers in liver injury and tumor progression could improve cancer diagnosis, prognosis, and management. To study this, we performed next generation sequencing of the livers of Sprague-Dawley rats treated with TAA at three doses (4.5, 15 and 45 mg/kg) and four time points (3-, 7-, 14- and 28-days). Overall, 330 unique differentially expressed miRNAs (DEMs) were identified in the entire TAA-treatment course. Of these, 129 DEMs were found significantly enriched for the "liver cancer" annotation. These results were further complemented by pathway analysis (Molecular Mechanisms of Cancer, p53-, TGF-ß-, MAPK- and Wnt-signaling). Two miRNAs (rno-miR-34a-5p and rno-miR-455-3p) out of 48 overlapping DEMs were identified to be early and sensitive biomarkers for TAA-induced hepatocarcinogenicity. We have shown significant regulatory associations between DEMs and TAA-induced liver carcinogenesis at an earlier stage than histopathological features. Most importantly, miR-34a-5p is the most suitable early and sensitive biomarker for TAA-induced hepatocarcinogenesis due to its consistent elevation during the entire treatment course.

Plasma miR-208 as a useful biomarker for drug-induced cardiotoxicity in rats.

Cardiotoxicity is one of the major safety concerns in drug development. Therefore, detecting and monitoring cardiotoxicity throughout preclinical and clinical studies is important for pharmaceutical companies. The present study was conducted in order to explore a plasma miRNA biomarker for cardiotoxicity in rats. As organ specificity is an important factor for a biomarker, we analyzed the miRNA microarray dataset in 55 organs/tissues in normal male rats. Based on this analysis, 5 miRNAs consisting of miR-208 (heart-specific), miR-1, miR-133a, miR-133b (heart and skeletal muscle-specific) and miR-206 (skeletal muscle-specific) were selected. Next, we evaluated the usefulness of those 5 miRNAs as circulating biomarkers in rats administered with single-dose isoproterenol or doxorubicin. Plasma miR-208 was consistently increased through 24 h after dosing in rats administered with isoproterenol, whereas plasma concentrations of cardiac troponin (cTn) showed transient elevation. In contrast, the plasma levels of miR-1, miR-133a, miR-133a and miR-206 were elevated after treatment with doxorubicin, probably as a result of skeletal muscle toxicity. Additionally, the plasma miR-208 level was elevated even after repeat-dose administration (once daily for 7 days) of isoproterenol under which the pathological condition proceeded to the sub-chronic phase such as fibrosis. Thus, our data suggest that miR-208 is a promising plasma biomarker for cardiotoxicity in rats. Monitoring of plasma miR-208 levels in rats may lead to more accurate evaluation of cardiotoxicity in preclinical studies.

Diagnostic and predictive performance and standardized threshold of traditional biomarkers for drug-induced liver injury in rats.

Traditional biomarkers such as alanine and aspartate aminotransferase (ALT, AST) and total bilirubin (TBIL) have been widely used for detecting drug-induced liver injury (DILI). Although the Food and Drug Administration (FDA) proposed standardized thresholds for human as Hy's law, those for animals have not been determined, and predictability of these biomarkers for future onset of hepatic lesions remains unclear. In this study, we investigated these diagnostic and predictive performance of 10 traditional biomarkers for liver injury by receiver-operating characteristic (ROC) curve, using a free-access database where 142 hepatotoxic or non-hepatotoxic compounds were administrated to male rats (n=5253). Standardization of each biomarker value was achieved by calculating the ratio to control mean value, and the thresholds were determined under the condition of permitting 5% false positive. Of these 10 biomarkers, AST showed the best diagnostic performance. Furthermore, ALT and TBIL also showed high performance under the situation of hepatocellular necrosis and bile duct injury, respectively. Additionally, the availability of the diagnostic thresholds in difference testing facility was confirmed by the application of these thresholds to in-house prepared dataset. Meanwhile, incorrect diagnosis by the thresholds was also observed. Regarding prediction, all 10 biomarkers showed insufficient performance for future onset of hepatic lesions. In conclusion, the standardized diagnostic thresholds enable consistent evaluation of traditional biomarkers among different facilities, whereas it was suggested that novel biomarker is required for more accurate diagnosis and prediction of DILI.

Comprehensive analysis of DNA methylation and gene expression of rat liver in a 2-stage hepatocarcinogenesis model.

Recent studies have shown that epigenetic alterations correlate with carcinogenesis in various tissues. Identification of these alterations might help characterize the early stages of carcinogenesis. We comprehensively analyzed DNA methylation and gene expression in livers obtained from rats exposed to nitrosodiethylamine (DEN) followed by a promoter of hepatic carcinogenesis, phenobarbital (PB). The combination of DEN and PB induced marked increases in number and area of glutathione S-transferase-placental form (GST-P)-positive foci in the liver. In the liver of rats that received 30 mg/kg of DEN, pathway analysis revealed alterations of common genes in terms of gene expression and DNA methylation, and that these alterations were related to immune responses. Hierarchical clustering analysis of the expression of common genes from public data obtained through the Toxicogenomics Project-Genomics Assisted Toxicity Evaluation system (TG-GATEs) showed that carcinogenic compounds clustered together. MBD-seq and GeneChip analysis indicated that major histocompatibility complex class Ib gene RT1-CE5, which has an important role in antigen presentation, was hypomethylated around the promoter region and specifically induced in the livers of DEN-treated rats. Further, immunohistochemical analysis indicated that the co-localization of GST-P and protein homologous to RT1-CE5 was present at the foci of some regions. These results suggest that common genes were altered in terms of both DNA methylation and expression in livers, with preneoplastic foci indicating carcinogenic potential, and that immune responses are involved in early carcinogenesis. In conclusion, the present study identified a specific profile of DNA methylation and gene expression in livers with preneoplastic foci. Early epigenetic perturbations of immune responses might correlate with the early stages of hepatocarcinogenesis.

Combinatorial measurement of CDKN1A/p21 and KIF20A expression for discrimination of DNA damage-induced clastogenicity.

In vitro mammalian cytogenetic tests detect chromosomal aberrations and are used for testing the genotoxicity of compounds. This study aimed to identify a supportive genomic biomarker could minimize the risk of misjudgments and aid appropriate decision making in genotoxicity testing. Human lymphoblastoid TK6 cells were treated with each of six DNA damage-inducing genotoxins (clastogens) or two genotoxins that do not cause DNA damage. Cells were exposed to each compound for 4 h, and gene expression was comprehensively examined using Affymetrix U133A microarrays. Toxicogenomic analysis revealed characteristic alterations in the expression of genes included in cyclin-dependent kinase inhibitor 1A (CDKN1A/p21)-centered network. The majority of genes included in this network were upregulated on treatment with DNA damage-inducing clastogens. The network, however, also included kinesin family member 20A (KIF20A) downregulated by treatment with all the DNA damage-inducing clastogens. Downregulation of KIF20A expression was successfully confirmed using additional DNA damage-inducing clastogens. Our analysis also demonstrated that nucleic acid constituents falsely downregulated the expression of KIF20A, possibly via p16 activation, independently of the CDKN1A signaling pathway. Our results indicate the potential of KIF20A as a supportive biomarker for clastogenicity judgment and possible mechanisms involved in KIF20A downregulation in DNA damage and non-DNA damage signaling networks.

Detection of initiating potential of non-genotoxic carcinogens in a two-stage hepatocarcinogenesis study in rats.

We previously reported a toxicogenomics-based prediction model for hepatocarcinogens in which the expression patterns of signature genes following repeated doses of either genotoxic or non genotoxic compounds were similar. Based on the results of our prediction model, we hypothesized that repeated doses of non-genotoxic carcinogens might have initiating potential. Here, we conducted a two stage hepatocarcinogenesis study in rats exposed to the initiating agent nitrosodiethylamine (DEN), and hepatotoxic compounds thioacetamide (TAA), methapyrilene (MP) and acetaminophen (APAP) for 1-2 weeks followed by the liver tumor promoter phenobarbital (PB). The duration of initial treatment was determined based on positive results from our prediction model. Combined treatment of 3 or 30 mg/kg of genotoxic DEN and PB induced marked increases in altered hepatocellular foci and a DEN dose-dependent increase in the number and area of glutathione S-transferase-placental form (GST-P)-positive foci. A low number of altered hepatocellular foci were also observed in rats treated with TAA at a dose of 45 mg/kg.MP at a dose of 100 mg/kg induced a very low number of foci, but APAP did not. Hierarchical clustering analysis using gene expression data revealed that 2-week treatment with TAA at a dose of 30 mg/kg and MP at 45 mg/kg induced specific expression of DNA damage-related genes, similar to 1-week treatment with DEN at a dose of 30 mg/kg. These results suggest that TAA and MP induce DNA damage, which partially supports our hypothesis. Although this study does not indicate whether tumor growth in response to these compounds can be assessed in this model, our results suggest that cumulative treatment with non genotoxic TAA might have initiating potential in the liver.

Comparative gene and protein expression analyses of a panel of cytokines in acute and chronic drug-induced liver injury in rats.

Drug-induced liver injury (DILI) is a significant safety issue associated with medication use, and is the major cause of failures in drug development and withdrawal in post marketing. Cytokines are signaling molecules produced and secreted by immune cells and play crucial roles in the progression of DILI. Although there are numerous reports of cytokine changes in several DILI models, a comprehensive analysis of cytokine expression changes in rat liver injury induced by various compounds has, to the best of our knowledge, not been performed. In the past several years, we have built a public, free, large-scale toxicogenomics database, called Open TG-GATEs, containing microarray data and toxicity data of the liver of rats treated with various hepatotoxic compounds. In this study, we measured the protein expression levels of a panel of 24 cytokines in frozen liver of rats treated with a total of 20 compounds, obtained in the original study that formed the basis of the Open TG-GATEs database and analyzed protein expression profiles combined with mRNA expression profiles to investigate the correlation between mRNA and protein expression levels. As a result, we demonstrated significant correlations between mRNA and protein expression changes for interleukin (IL)-1ß, IL-1α, monocyte chemo-attractant protein (MCP)-1/CC-chemokine ligand (Ccl)2, vascular endothelial growth factor A (VEGF-A), and regulated upon activation normal T cell expressed and secreted (RANTES)/Ccl5 in several different types of DILI. We also demonstrated that IL-1ß protein and MCP-1/Ccl2 mRNA were commonly up-regulated in the liver of rats treated with different classes of hepatotoxicants and exhibited the highest accuracy in the detection of hepatotoxicity. The results also demonstrate that hepatic mRNA changes do not always correlate with protein changes of cytokines in the liver. This is the first study to provide a comprehensive analysis of mRNA-protein correlations of factors involved in various types of DILI, as well as additional insights into the importance of understanding complex cytokine expression changes in assessing DILI.

Identification of metabolomic biomarkers for drug-induced acute kidney injury in rats.

Nephrotoxicity is a common side effect observed during both nonclinical and clinical drug development investigations. The present study aimed to identify metabolomic biomarkers that could provide early and sensitive indication of nephrotoxicity in rats. Metabolomic analyses were performed using capillary electrophoresis-time-of-flight mass spectrometry on rat plasma collected at 9 and 24 h after a single dose of 2-bromoethylamine or n-phenylanthranilic acid and at 24 h after 7 days of repeated doses of gentamicin, cyclosporine A or cisplatin. Among a total of 169 metabolites identified, 3-methylhistidine (3-MH), 3-indoxyl sulfate (3-IS) and guanidoacetate (GAA) were selected as candidate biomarkers. The biological significance and reproducibility of the observed changes were monitored over time in acute nephrotoxicity model rats treated with a single dose of cisplatin, with the glomerular filtration rate monitored by determination of creatinine clearance. Increased plasma levels of 3-MH and 3-IS were related to a decline in glomerular filtration due to a renal failure. In contrast, the decrease in plasma GAA, which is synthesized from arginine and glycine in the kidneys, was considered to reflect decreased production due to renal malfunction. Although definitive validation studies are required to confirm their usefulness and reliability, 3-MH, 3-IS and GAA may prove to be valuable plasma biomarkers for monitoring nephrotoxicity in rats.

miRNA expression atlas in male rat.

MicroRNAs (miRNAs) are small (~22 nucleotide) noncoding RNAs that play pivotal roles in regulation of gene expression. The value of miRNAs as circulating biomarkers is now broadly recognized; such tissue-specific biomarkers can be used to monitor tissue injury and several pathophysiological conditions in organs. In addition, miRNA profiles of normal organs and tissues are important for obtaining a better understanding of the source of modulated miRNAs in blood and how those modulations reflect various physiological and toxicological conditions. This work was aimed at creating an miRNA atlas in rats, as part of a collaborative effort with the Toxicogenomics Informatics Project in Japan (TGP2). We analyzed genome-wide miRNA profiles of 55 different organs and tissues obtained from normal male rats using miRNA arrays. The work presented herein represents a comprehensive dataset derived from normal samples profiled in a single study. Here we present the whole dataset with miRNA profiles of multiple organs, as well as precise information on experimental procedures and organ-specific miRNAs identified in this dataset.

Genomic biomarkers for cardiotoxicity in rats as a sensitive tool in preclinical studies.

The development of safer drugs is a high priority for pharmaceutical companies. Among the various toxicities caused by drugs, cardiotoxicity is an important issue because of its lethality. In addition, cardiovascular toxicity leads to the attrition of many drug candidates in both preclinical and clinical phases. Although histopathological and blood chemistry examinations are the current gold standards for detecting cardiotoxicity in preclinical studies, the large number of withdrawals from clinical studies owing to safety problems indicate that a more sensitive tool is required. We recently identified 32 genes that were candidate genomic biomarkers for cardiotoxicity in rats. Based on their functions, the present study focused on 8 of these 32 genes (Spp1, Fhl1, Timp1, Serpine1, Bcat1, Lmcd1, Rnd1 and Tgfb2). Diagnostic accuracy for the genes was determined by a receiver-operating characteristic (ROC) analysis using more cardiotoxic and non-cardiotoxic compounds. In addition, an optimized support vector machine (SVM) model that was composed of Spp1 and Timp1 was newly constructed. This new multi-gene model exhibited a much higher diagnostic accuracy than that observed for plasma cardiac troponin I (cTnI), which is one of the most useful plasma biomarkers for cardiotoxicity detection. Furthermore, we determined that this multi-gene model could predict potential cardiotoxicity in rats in the absence of any cardiac histopathological lesions or elevations of plasma cTnI. Overall, this multi-gene model exhibited advantages over classic tools commonly used for cardiotoxicity evaluations in rats. Our current results suggest that application of the model could potentially lead to the production of safer drugs.

Toxicogenomics discrimination of potential hepatocarcinogenicity of non-genotoxic compounds in rat liver.

Long-term carcinogenicity testing of a compound is exceedingly time-consuming and costly, and requires many test animals, whereas the Ames test, which is based on the assumption that any substance that is mutagenic may also exert carcinogenic potential, is useful as a short-term screening assay but has major drawbacks. Although, in fact, 90% of compounds that give a positive Ames test cause cancer in laboratory animals, a good proportion of compounds that give a negative Ames test are also carcinogens; that is, there is no good correlation between carcinogenicity and negative Ames test results. As an alternative to these two approaches, we have tried applying toxicogenomics to predict the carcinogenicity of a compound from the gene expression profile induced in vivo. To establish our model, male Sprague-Dawley rats were orally administered test compounds (12 hepatocarcinogens and 26 non-hepatocarcinogens) for 28 days. Analysis of liver gene expression data by Support Vector Machines (SVM) dividing compounds into 'for training' and 'for test' (20 cases assigned randomly) allowed a set of marker genes to be tested for prediction of hepatocarcinogenicity. The developed prediction model was then validated with reference to the concordance rate with training data and test data, and a good performance was obtained. We will have new gene expression data and continue the validation of our model.

Toxicogenomic biomarkers for renal papillary injury in rats.

Renal papillary injury is a common side effect observed during nonclinical and clinical investigations in drug development. The present study aimed to identify genomic biomarkers for early and sensitive detection of renal papillary injury in rats. We hypothesized that previously identified genomic biomarkers for tubular injury might be applicable for the sensitive detection of papillary injury in rats. We selected 18 genes as candidate biomarkers for papillary injury based on previously published studies and analyzed their expression profiles by RT-PCR in each kidney region, namely the cortex, cortico-medullary junction, and papilla in various nephrotoxicity models. Comparative analysis of gene expression profiles revealed that some genes were commonly upregulated or downregulated in the renal papilla, reflecting papillary injuries induced by 2-bromoethylamine hydrobromide, phenylbutazone, or n-phenylanthranilic acid. By applying receiver operator characteristics analysis, six candidate biomarkers were identified and their usefulness was confirmed by using an independent data set. The three top-ranked genes, Timp1, Igf1, and Lamc2, exhibited the best prediction performance in an external data set with area under the curve (AUC) values of greater than 0.91. An optimized support vector machine model consisting of three genes achieved the highest AUC value of 0.99. In conclusion, even though definitive validation studies are required for the establishment of their usefulness and reliability, these identified genes may prove to be the most promising candidate genomic biomarkers of renal papillary injury in rats.

Underestimation of urinary biomarker-to-creatinine ratio resulting from age-related gain in muscle mass in rats.

Recent efforts have been made to identify useful urinary biomarkers of nephrotoxicity. Furthermore, the application of urine to the other toxicities as new biomarker source has been recently expanded. Meanwhile, correction of urinary biomarker concentrations according to fluctuations in urine flow rate is required for adequate interpretation of the alteration. The urinary biomarker-to-creatinine ratio (UBCR) is widely used because of the convenience, while the urinary biomarker-excretion rate is regarded as the gold standard corrective method. Because creatinine is a catabolite in energy production in muscles, we hypothesized that altered muscle mass could affect creatinine kinetics, ultimately affecting UBCR. However, no study has examined this hypothesis. In this study, we examined the influence of muscle mass gain on UBCR, using male Sprague-Dawley rats during the growth phase, 6-12-week old. Both plasma creatinine and excretion of urinary creatinine (Ucr excretion) showed increases with muscle mass gain in rats, in which the alterations of UBCR were lowered. The renal mRNA level of the organic cation transporter-2 (Oct2), a creatinine transporter, showed an age-related increase, whereas the mRNA level of multidrug and toxin extrusions-1 (Mate1) remained constant. Multiple regression analysis showed that the increase in creatinine clearance highly contributed to the age-related increase in Ucr excretion compared to the mRNA levels of Oct2 and Mate1. This suggested that the age-related increase in Ucr excretion may be attributable to the increased transglomerular passage of creatinine. In conclusion, the results suggest that muscle mass gain can affect creatinine kinetics, leading to underestimation of UBCR. Therefore, it is important to understand the characteristics of the corrective method when using urinary biomarker, the failure of which can result in an incorrect diagnosis.

Toxicogenomic multigene biomarker for predicting the future onset of proximal tubular injury in rats.

Drug-induced renal tubular injury is a major concern in the preclinical safety evaluation of drug candidates. Toxicogenomics is now a generally accepted tool for identifying chemicals with potential safety problems. The specific aim of the present study was to develop a model for use in predicting the future onset of drug-induced proximal tubular injury following repeated dosing with various nephrotoxicants. In total, 41 nephrotoxic and nonnephrotoxic compounds were used for the present analysis. Male Sprague-Dawley rats were dosed orally or intravenously once daily. Animals were exposed to three different doses (low, middle, and high) of each compound, and kidney tissue was collected at 3, 6, 9, and 24 h after single dosing, and on days 4, 8, 15, and 29 after repeated dosing. Gene expression profiles were generated from kidney total RNA using Affymetrix DNA microarrays. Filter-type gene selection and linear classification algorithms were employed to discriminate future onset of proximal tubular injury. We identified genomic biomarkers for use in future onset prediction using the gene expression profiles determined on day 1, when most of the nephrotoxicants had yet to produce detectable histopathological changes. The model was evaluated using a five-fold cross validation, and achieved a sensitivity of 93% and selectivity of 90% with 19 probes. We also found that the prediction accuracy of the optimized model was substantially higher than that produced by any of the single genomic biomarkers or histopathology. The genes included in our model were primarily involved in DNA replication, cell cycle control, apoptosis, and responses to oxidative stress and chemical stimuli. In summary, our toxicogenomic model is particularly useful for predicting the future onset of proximal tubular injury.

Prediction model of potential hepatocarcinogenicity of rat hepatocarcinogens using a large-scale toxicogenomics database.

The present study was performed to develop a robust gene-based prediction model for early assessment of potential hepatocarcinogenicity of chemicals in rats by using our toxicogenomics database, TG-GATEs (Genomics-Assisted Toxicity Evaluation System developed by the Toxicogenomics Project in Japan). The positive training set consisted of high- or middle-dose groups that received 6 different non-genotoxic hepatocarcinogens during a 28-day period. The negative training set consisted of high- or middle-dose groups of 54 non-carcinogens. Support vector machine combined with wrapper-type gene selection algorithms was used for modeling. Consequently, our best classifier yielded prediction accuracies for hepatocarcinogenicity of 99% sensitivity and 97% specificity in the training data set, and false positive prediction was almost completely eliminated. Pathway analysis of feature genes revealed that the mitogen-activated protein kinase p38- and phosphatidylinositol-3-kinase-centered interactome and the v-myc myelocytomatosis viral oncogene homolog-centered interactome were the 2 most significant networks. The usefulness and robustness of our predictor were further confirmed in an independent validation data set obtained from the public database. Interestingly, similar positive predictions were obtained in several genotoxic hepatocarcinogens as well as non-genotoxic hepatocarcinogens. These results indicate that the expression profiles of our newly selected candidate biomarker genes might be common characteristics in the early stage of carcinogenesis for both genotoxic and non-genotoxic carcinogens in the rat liver. Our toxicogenomic model might be useful for the prospective screening of hepatocarcinogenicity of compounds and prioritization of compounds for carcinogenicity testing.